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We present the computation, in lattice QCD, of the renormalization constants and 
mixing coefficients of operators that measure the first two moments of DIS Struc- 
ture Functions. These calculations have been performed using the Sheikholeslami- 
Wohlert O(a) improved "clover" action, which is known to reduce the systematic 
error associated with the finiteness of the lattice spacing a. Due to the com- 
plexities of the computations, we have developed, using the computer languages 
Schoonschip and Form, general codes that are able to automatically carry out all 
the analytic lattice manipulations. 



1 Introduction 

The computation on the lattice of renormalization constants is a necessary 
ingredient for the connection of lattice operators and matrix elements to their 
continuum counterparts and the extraction of physical quantities from the 
numbers obtained in Monte Carlo simulations. We have computed, in lattice 
1-loop perturbation theory using the Sheikholeslami-Wohlert O(a) improved 
action, the renormalization constants and mixing coefficients of the quark an 
gluon operators of rank twoEJu and of the quark operators of rank three □ 
that measure n th£ first two moments of DIS Structure Functions. The use of 
improvement au allows us to reduce the systematic error associated with the 
finiteness of the lattice spacing a. In the particular formulation we have used 
this error is lowered, for on-shcll quantities, from 0(a) to 0(a/loga). 

2 Moments of Structure Functions 

The hadronic tensor W MIy , from which the Structure Functions can be defined 
in a well-known way, is written in terms of the hadronic currents as 

= 77" fdtx e^ x < p\J^{x)J v {Q)\p > . (1) 
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By means of a Wilson OPE expansion near the light-like region of the kind 
A(x)B(0) ~ J2n i c N,i{x 2 ) x^ 1 ■ ■ ■ x tJ,N Of^' i ) flN (0) , it can be expressed in terms 

of a set of symmetric and traceless operators Off^.'j)^ with vanishing vacuum 
expectation values, of which the dominant ones have twist two. The matrix 
elements of these operators have the general form 



< p\0}£).. m \p >= Ajv(^)p w • • ■p f _ lN + trace terms, 
and are related to the moments of the Structure Functions by the formula 



(2) 



dxx N - 1 ^ k (q 2 ,x) = C N (q 2 /fJ, 2 )A N (fi), 



(3) 



where T\ — 1F\ , T2 = -F2 /x and JF3 = F3 . The coefficients Cn are known 
from continuum P.T., and thus we can extract a given moment < x N > once 
we know the corresponding matrix element < p\0^1.^ N \p >. Such matrix 
elements contain long distance (non-perturbative) physics, thus the only viable 
way to compute the moments of Structure Functions is with the use of lattice 
methods. 

We have considered in our calculations the unpolarized Structure Func- 
tions, and in particular we have computed the operators 



O q 



1 



< x > Q 



< x >, 



(4) 
(5) 



O q 



< X* > 



(6) 



3 Improved Lattice QCD 



Lattice QCD allows the evaluation from first principles of the hadronic matrix 
elements needed for the computation of the moments of the Structure Func- 
tions. Once discretization is introduced, the quark fields, ip n , exist on the 
lattice sites, and the gauge fields, U„^ = e t9 " at A ™,», exist as links between 
these sites. The Wilson actionem 



- i> n ( m f + ) ip n 



9l 
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2 



is widely used as discretization of the (Euclidean) QCD action. The corre- 
sponding regularization scheme is gauge- invariant, but the terms proportional 
to r (0 < r < 1), introduced to get rid of lattice spurious fermions, break chiral 
symmetry even in the case of a vanishing quark mass m/. 

A matrix element like < plO^)..^ \p > (i.e. a given moment of the Struc- 
ture Functions) can be determined from the computation of two- and three- 
point correlation functions. However, its determination from Monte Carlo sim- 
ulations is affected by both statistical and systematic errors. The statistical 
errors come from the finite number of configurations used, while the systematic 
errors are of various nature: finite lattice spacing a, finite volume V, quenched 
approximation (that is, dropping the contribution of the internal quark loops), 
and from the necessary extrapolation to recover chiral symmetry. 

We are interested here in the systematic error associated to the finiteness 
of the lattice spacing, and to reduce this error improvement procedures have 
been proposedtT LJ The formulation that we use consists of adding an irrelevant 
operator to the standard Wilson action in such a way to cancel, in on-shell 
matrix elementsQ all terms that in the continuum limit jare effectively of order 
"a" . This "clover-leaf" Shcikholeslami- Wohlert term is □ 

ASf = -iffoa 4 ^2 ^ ' t Pn cr ^ F n,^1pn, (8) 

where F n ^ v is the average of the four plaquettes lying in the plane jiv, stem- 
ming from the point n: 

Fn,nu = T ^ Pn.^v = „. 2 {Un,fiu - (9) 

with U n ^,j = Un^Un+^uUn+v JJn,w This new term in the action means that 
we have to add to the Wilson quark-quark-gluon interaction vertex, 



(vt c = - go (t A ) bc 



a{k + k') p a{k + k')„ 
r sin h ij p cos 



(10) 



the improved quark-quark-gluon interaction vertex 

(V 1 )"; = -9^{t A ) bc cos °( fc - fc, )p J- g p\ sina(fc - k%. (11) 

A 

The fermion propagator is not modified by this improved action, and neither 
is the gluon propagator as the first corrections to the pure gauge term of the 
Wilson action are already of order a 2 . 
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As well as adding the term (g), in the calculation of a fermionic Green 
function the fermion fields undergo the rotation 



'0 



if) ► 1p 



ar 

1 + T 



(12) 



This rotation combined with the use of the Sheikholeslami-Wohlert action has 
been shown to remove, from on-shell hadronic matrix .elements, all terms that 
in the continuum limit are effectively of order autB Using this recipe, the 
systematic error related to the lattice discretization drops from order a to 
order a/ log a: 



P 



O, 



P 



Monte Carlo 







[(■> 


d 



p 



phys 



+ 0(a/ log a) 



(13) 



The magnitude of the order a terms is about 20-30 %, while the magni- 
tude of the order a/ log a terms is about 5-10 %. Therefore, with the use of 
Sheikholeslami-Wohlert improvement one can achieve a remarkable decrease 
of the systematic error coming from the finiteness of the lattice spacing. 



4 Renormalization constants 

The renormalization constants connect the bare lattice operators, O(a), to 
finite operators, 0(/i), renormalized at a scale ^i: 

O l (fi) = Z lk (fia)O k (a). (14) 

These constants are fixed in perturbation theory by the same renormalization 
conditions used in the continuum. In the flavor Singlet case there is a mix- 
ing between quark and gluon operators of the same rank that have the same 
conserved quantum numbers. We can then write: 

0« = Z qq Oi + z qg o* 

ds = z gq oi + z gg os, v ' 

and in this case all elements of the mixing matrix 

< q\CP\q > < g,<j\Oi\g,cr> \ 

< q\03\q > < g,a\03\g,n> ) K ' 

have to be computed. 

To this mixing (already present in continuum QCD) the lattice formula- 
tion adds additional mixings, induced by the breaking of (Euclidean) Lorentz 
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invarianceEil In some cases it is possible, by a careful choice of the Lorentz 
indices, to choose operators that are multiplicatively renormalizable on the 
lattice, but the higher the moment the more complicated the mixing pattern 
isB13 The breaking of the Lorentz invariance has also bound us to develop 
special computer routines to correctly perform the Dirac algebra on the lat- 
tice. They are a necessary ingredient in our general codes that are able to 
automatically carry out all the stages of the algebraic manipulations from the 
elementary building blocks of each Feynman diagram. 



5 Some results 



Some sinmlations have been performed in the past with the unimproved Wilson 
actionEMj Within errors, the results of these simulations are consistent with 
experiment. However, the values of the renormalization constants change non- 
trivially when they are computed in the improved theory. As an example we 
give here a selection of the results for quark operators. At j3 = 6/gg = 6.0 (the 
general dependence is Z = 1 + const//?), and for r — 1, we have: 

0? 12} = 1.027 O q Wilson 
(O q {n} y= 1.134 {O q {12} f Improved 

0\ 123} = 1-160 Of 123} Wilson 
{O q {im ) 1 = 1.252 (O q {123} ) J Improved 

Odis = I [1-184 A + 1.156 Ob] Wilson 
(Odis) 1 = I t 1 - 331 (OaY + 1-187 (Ob) 1 ] Improved, 



where Op IS = 0' 411 j — 5(^(422} + ^{433}) (which can be writtten as one-third 
of the sum of the non-symmetric operators Oa = 0\ xl — \ {0\ 2 2 + ^433) an( l 
Ob = 141 + O q u - \{0\ i2 + 0| 24 + 0| 43 + 0| 34 )) is not multiplicatively 
renormalizable on the lattice. We see there is a great difference between the 
Wilson and the improved results; in particular, the improved renormalization 
constants are somewhat higher. For this reason, improved simulations using 
the new values of the renormalization constants will give better insight into 
the agreement with the experimental values of the moments. 

Finally, we want to mention ithat some results are now available also for 
the polarized Structure Functionspl although so far limited to the Wilson case. 
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